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ABSTRACT 

We present the results of a multiline survey at mm-wavelengths of the Cepheus 
A star forming region. Four main flows have been identified: three pointing in the 
SW, NE, and SE directions and accelerating high density CS clumps. The fourth out- 
flow, revealed by high-sensitivity HDO observations, is pointing towards South and is 
associated with conditions particularly favourable to a chemical enrichment. At the 
CepA-East position the emissions due to the ambient clump and to the outflows coex- 
ist and different molecules exhibit different spectral behaviours. Some species (C 13 CH, 
C3H2, CH2CO, CH3C2H, HC 18 + ) exhibit relatively narrow lines at ambient veloc- 
ities (ambient peak). Other molecules (CO, CS, H 2 S, SiO, SO, SO2) show extended 
wings tracing the whole range of the outflow velocities. Finally, OCS, H 2 CS, HDO, and 
CH3OH are associated with wings and, in addition, show wings and in addition reveal 
a bright high velocity redshifted spectral peak (outflow peak) which can be used to 
investigate the southern outflows. At ambient velocities the gas is dense (> 10 5 cm~ 3 ) 
and different components at distinct temperatures coexist, ranging from the relatively 
low kinetic temperatures (< 50 K) measured with H2S, CH3OH, H2CS, and CH3C2H, 
to definitely higher temperature conditions, ~ 100-200 K, obtained from the SiO, SO, 
and SO2 spectra. For the outflow peak we derive densities between ~ 10 4 cm~ 3 to 
~ 10 7 cm~ 3 and high temperatures, ~ 100-200 K, indicating regions compressed and 
heated by shocks. 

The analysis of the line profiles shows that the SiO molecule dominates at the 
highest velocities and at the highest excitation conditions, confirming its close associ- 
ation with shocks. H2S, SO2, and SO preferentially trace more quiescent regions than 
SiO, and in particular a lack of bright H2S emission at the highest velocities is found. 
OCS and H2CS emit at quite high velocities, where the abundances of three shock 
tracers like SiO, CH3OH, and HDO arc higher. These results may indicate that H2S is 
not the only major sulphur carrier in the grain mantles, and that OCS and H2CS may 
probably play an important role on the grains; or that alternatively they rapidly form 
once the mantle is evaporated after the passage of a shock. Finally, the outflow peak 
emission has been compared with recent timc-dcpcndcnt sulphur chemistry models: 
the results indicate that, if associated with accurate measurements of the physical con- 
ditions, the CH3OH/H2CS column density ratio can be used as an effective chemical 
clock to date the age of shocked gas. 

Key words: ISM: clouds - ISM: individual objects: CepA - ISM: jets and outflows 
- ISM: molecules - Radio lines: ISM 



1 INTRODUCTION 

The mass loss from Young Stellar Objects (YSOs) produces 
high velocity flows which strike the ambient medium driv- 



ing shocks. Once the temperature has increased, at least 
to a few thousands degrees depending on the shock type, 
the energy barriers between neutral molecules can be over- 
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come and the chemistry of certain species, such as the sul- 
phuretted species, is altered significantly (see e.g. Pineau des 
Forets et al. 1993, van Dishoeck & Blake 1998, and reference 
therein). In addition, grains are affected by shocks, with the 
consequent injection of molecular and atomic species in the 
gas phase: again, this leads to an enhancement of the abun- 
dances of several species, including S-bearing molecules. The 
scenario proposed by most models is that HbS is the main 
reservoir of sulphur on grain mantles: once in the gas phase, 
H2S is used for a fast production of SO and SO2 (e.g. Pineau 
des Forets et al. 1993, Charnley 1997). However, the lack of 
H2S features in the ISO spectra (Gibb et al. 2000, Boogert 
et al. 2000) which set upper limits on the iced H2S abun- 
dance around protostars and the detection of OCS on grains 
(Palumbo et al. 1997) suggest that the latter may be an im- 
portant sulphur carrier in the ices. This seems supported 
by the observations in the envelopes of massive young stars 
recently performed by van der Tak et al. (2003), which indi- 
cate for OCS higher excitation temperatures than for H2S. 
An alternative hypothesis is that the sulphur released from 
the dust mantles is mainly in atomic form (Wakelam et al. 
2004). In any case, once the gas phase has been enriched by 
the passage of a shock, other S-bearing species such as H2CS 
and HCS + are expected to significantly increase their abun- 
dances as a consequence of the sulphur injection. Therefore, 
estimates of abundance ratios such as SO2/H2S, H2S/OCS, 
and H2CS/OCS may provide us with chemical clocks to 
study the evolutionary stages of molecular outflows. In fact, 
the use of SO2/H2S and H2S/OCS ratios has already led 
crude age estimates of the outflows located in CB3 (Codella 
& Bachiller 1999) and LI 157 (Bachiller et al. 2001), encour- 
aging further studies. 

The high temperature scenario can be applied also to 
hot cores around protostars, which are characterised by high 
temperatures (> 100 K), high densities (> 10 6 cm -3 ), and 
a rich molecular inventory. The SO2/H2S ratio has been 
used by Hatchell et al. (1998) for massive hot cores: the 
inferred ages are in agreement with the dynamical times 
estimated from the associated outflows. Moreover, Buckle 
& Fuller (2003) show that the chemical evolution of sul- 
phuretted species is a potential probe of timescales also in 
low-mass star forming regions. In summary, the effect of 
high-temperature chemistry on the composition of the gas 
hosting the star forming process can be used as a tool to 
investigate the evolution of protostars. 

Physical conditions as well as time evolution affect the 
sulphur chemistry in high-temperatures environment. For 
instance, recently Wakelam et al. (2004) modeled the sul- 
phur chemistry around hot core-like environments and found 
that the obtained abundances depend not only on age but 
also on the excitation conditions of the gas. Thus, in order 
to use the sulphur abundance ratios as chemical clocks, care 
should be taken to first constrain the gas conditions. In the 
case of molecular outflows, it is reasonable to expect different 
physical and/or chemical conditions at the different veloci- 
ties and thus a study of the excitation as a function of the 
velocity is needed. The existing observations of sulphuretted 
molecules have been carried out to date in a unsystematic 
way, and, in particular, since the wing profiles are weak, 
the chemical composition of the gas at high-velocity has 
been poorly investigated. With this in mind, we carried out 
line observations of S-bearing species at mm- wavelengths to- 



wards a star forming region with well defined high-velocity 
components. Cepheus A (CepA) well represents such target: 
in particular, the East component harbours an OB3 stellar 
association (Goetz et al. 1998, and references therein) driv- 
ing multiple outflows (e.g. Narayanan & Walker 1996, Bergin 
et al. 1997). Bergin et al. (1997) have studied the chemi- 
cal properties of Cepheus A showing that a large number of 
molecules can be easily detected, but their observations only 
trace the quiescent gas, this is likely due to the limitation 
in spatial and spectral resolutions, which does not allow the 
detection of the high velocity wings. 

In this paper, we report the results of a deep millimeter 
survey of molecular lines in CepA. In particular, we compare 
line observations of SO, HCS+, H 2 CS, OCS, H 2 S, and SO2 
with profiles due to standard tracers of shocked gas, high 
excitation conditions, and of ambient emission. The main 
aim is to study the variation of the abundance ratios and 
the excitation conditions of S-bearing species along the line 
profiles. We find that sulphuretted molecules may be good 
chemical clocks to date molecular outflows and their driving 
protostellar sources. 



2 OBSERVATIONS AND RESULTS 

The observations were performed with the IRAM 30-m tele- 
scope at Pico Veleta (Granada, Spain) in June 2001, Septem- 
ber 2002, August 2003, and June 2004. The observed molec- 
ular species, the transitions, their upper level energies, the 
Einstein A coefficients and rest frequencies and some observ- 
ing parameters, such as the HPBW and the typical system 
temperature (T sys ) are summarised in Table 1. The integra- 
tion time (ON+OFF source) ranged from about 1 to ~18 
hours, depending on the intensities of the observed lines. 
The main beam efficiency varies from about 0.8 (at 81 GHz) 
to 0.5 (at 262 GHz). The observations were made by posi- 
tion switching. The pointing was checked about every hour 
by observing nearby planets or continuum sources and it 
was found to be accurate to within 4". As spectrometers, 
an autocorrelator (AC) split into different parts (up to six) 
was used to allow simultaneous observations of four different 
transitions. Also a 1 MHz filter bank, split into four parts of 
256 channels, was simultaneously used. The velocity resolu- 
tions provided by both backends, AC and 1 MHz, are shown 
in Table 1. When considered convenient, the AC spectra 
were smoothed to a lower velocity resolution (up to ~ 1 km 
s _1 ). The spectra were calibrated with the standard chopper 
wheel method and are reported here in units of main-beam 
brightness temperature (Tmb). 

The multiline molecular survey has been performed to- 
wards CepA-East and, in particular, at the coordinates of 
the HW2 object, which is one of the YSOs thought to drive 
the molecular outflows: a 2 ooo = 22 h 56 m 1719, foooo = +62° 
01' 49'.'7. In addition, we present maps of CepA-East in CS, 
a well known tracer of high density clumps, and in HDO, 
which traces hot gas chemistry as well as shocked material 
(e.g. van Dishoeck & Blake 1998). In particular, the heavy 
water emission will be here used as alternative to H2O, 
whose observations are prevented from the ground due to 
very strong atmospheric absorption. Water is expected to 
be one of the most important coolant in non-dissociative 
shocks (e.g. Kaufman & Neufeld 1996), since its gas-phase 
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Table 1. List of molecular species, transitions and observing parameters 



Transition 


u o 


Eu A u i F int 


HPBW T sys 


dv(AC) 


(fo(lMHz) 




(MHz) 


(K) (s- 1 ) (K km/s) 


(") (K) 


(km/s) 


(km/s) 


Selected transitions 



HDO(J k _k + = lio-lii) 


80578.30 


47 


1.3 10 


-6 


0.58(0.05) 


31 


150 


0.15 


3.72 


OCS(J = 7-6) 


85139.12 


16 


1.7 10 


-6 


0.36(0.02) 


29 


110 


0.14 




HCS+(J = 2-1) 


85347.88 


6 


1.1 10 


-5 


0.38(0.04) 


29 


110 


0.14 


3.51 


CH 3 C 2 H(Jk = 5 K -4 K ) 


85457.30° 


12° 


2.0 10" 


6 a 


3.11(0.13) ft 


29 


150 


0.27 


3.51 


C 34 S(J = 2-1) 


96412.98 


7 


1.6 10 


-5 


2.16(0.18) 


26 


220 


0.24 


3.11 


CS(J = 2-1) 


97980.97 


7 


1.7 10 


-5 


28.03(0.14) 


25 


180 


0.24 


3.06 


SO(J K = 3 2 -2i) 


99299.88 


4 


1.1 10 


-5 


14.42(0.03) 


24 


150 


0.24 


3.02 


H 2 CS(Jk_k + — 3i3— 212) 


101477.75 


23 


1.3 10 


-5 


1.01(0.03) 


24 


187 


0.12 


2.95 


34 S0 2 (Jk_k + = 3i 3 -2 02 ) 


102031.91 


8 


9.5 10 


-6 


0.11(0.01) 


24 


130 


0.11 


2.94 


CH 3 C 2 H(J K = 6 K -5 K ) 


102547.98" 


17° 


3.4 10" 


6a 


4.16(0.14) ft 


24 


200 


0.23 


2.90 


SO(J K = 4 3 -3 2 ) 


138178.64 


9 


3.1 10 


-5 


22.76(0.19) 


17 


370 


0.17 


2.17 


CH 3 OH(J K = 3 K -2 K ) 


145103.23° 


14° 


1.2 10" 


5a 


5.64(0.15)° 


17 


280 


0.08 


2.07 


H 2 34 S(J k _k + = lio-loi) 


167910.52 


28 


2.6 10 


-5 


1.11(0.23) 


14 


480 


0.07 


1.79 


H 2 CS(Jk_k + = 6i6-5i 5 ) 


202923.55 


47 


1.2 10 


-4 


2.15(0.08) 


12 


460 


0.06 


1.48 


CH 3 C 2 H(J K = 12 k -11k) 


205080.73° 


64° 


2.8 10" 


5a 


6.71(0.36) 6 


12 


600 


0.11 


1.46 


HCS+(J = 5-4) 


213360.64 


31 


2.0 10 


-4 


0.59(0.08) 


12 


330 


0.04 


1.41 


H 2 34 S(J k _k + = 2 20 -2n) 


213376.92 


84 


1.7 10 


-5 


<0.07 c 


12 


330 


0.04 


1.41 


SiO(J = 5-4) 


217104.94 


31 


5.2 10 


-4 


8.71(0.24) 


11 


450 


0.43 


1.38 


C ls O(J = 2-1) 


219560.33 


16 


6.2 10 


-7 


47.94(0.25) 


11 


580 


0.43 


1.36 


SO(J K = 6 5 -5 4 ) 


219949.39 


24 


1.3 10 


-4 


40.74(0.30) 


11 


580 


0.11 


1.36 


CO(J = 2-1) 


230537.98 


17 


6.9 10 


-7 




10 


1100 


0.10 


1.30 


C 34 S(J = 5-4) 


241016.17 


34 


2.8 10 


-4 


1.75(0.40) 


10 


990 


0.10 


1.24 


HDO(J K k + = 2n-2 12 ) 


241561.53 


95 


1.2 10 


-5 


3.42(0.19) 


10 


68 


0.05 


1.24 


CH 3 OH(J K = 5 K -4 K ) 


241791.44° 


35° 


5.8 10" 


5a 


8.32(0.07)° 


10 


678 


0.05 


1.24 


H 2 CS(Jk_k + = 7i6-6is) 


244047.75 


60 


2.1 10 


-4 


1.36(0.06) 


10 


890 


0.05 


1.23 


CS(J = 5-4) 


244935.61 


34 


3.0 10 


-4 


47.33(0.53) 


10 


1100 


0.10 


1.22 


SO(J K = 7 6 -6 5 ) 


261843.72 


35 


2.2 10 


-4 


52.30(0.58) 


9 


1300 


0.09 


1.15 



Serendipity detections 



HC ls O+(J = 1-0) 


85162.21 


4 


3.6 10" 


-5 


0.65(0.02) 


29 


110 


0.14 




C 13 CH(J K _ K+ = lii-Oii) 


85307.69 


4 


1.3 10" 


-6 


0.18(0.03) 


29 


110 




3.51 


C 3 H 2 (Jk_k + = 2i 2 -loi) 


85338.91 


6 


2.3 10" 


-5 


1.18(0.05) 


29 


110 


0.14 


3.51 


CH 2 CO(J k _k + = 5i4-4i 3 ) 


101981.43 


23 


1.1 10" 


-5 


0.14(0.01) 


24 


130 


0.11 


2.94 


CH 3 OH(J K = 10_2-10i E) 


102122.70 


154 


1.7 10" 


-7 


0.08(0.01) 


24 


130 




2.94 


CH 3 OH(J K = 9i-9 E) 


167931.13 


126 


2.3 10" 


-5 


2.33(0.25) 


14 


480 


0.07 


1.79 


CH 3 OH(J K = 13 6 -14 5 E) 


213377.52 


390 


1.1 10" 


-5 


0.43(0.03) 


12 


330 


0.04 


1.41 


CH 3 OH(J K = li-Oo E) 


213427.12 


13 


3.4 10" 


-5 


2.20(0.06) 


12 


330 


0.04 


1.41 


34 S0 2 (J k _k+ = 16ii5-15 2 i 4 ) 


241509.05 


131 


8.3 10" 


-5 


0.78(0.14) 


10 


670 




1.24 


HNCO(J k _k + = llon-lOoio) 


241774.09 


58 


2.0 10" 


-4 


2.81(0.20) 


10 


678 


0.05 


1.24 


34 S0 2 (J k _k + = I8117-I8018) 


243935.88 


163 


7.0 10" 


-5 


0.61(0.73) 


10 


890 




1.23 



° For CH 3 OH, it refers to the Jo-J-lo A+ line, while for CH 3 C 2 H it refers to the K=0 component. b For CH 3 C 2 H the integrated flux 
refers to the K=0,1 pattern. c It refers to the 3cr noise level; d The CO(2-l) observations have been performed in wobbler mode loosing 
the information at ambient velocities in order to investigate the high- velocity wings. 



abundance is considerably enhanced both via high temper- 
ature (larger than 200-300 K) reactions and by sputtering. 
Given the characteristics of CepA-East, it will thus be pos- 
sible to examine in details how the dense medium is affected 
by the occurrence of YSOs 

Table 1 reports the selected transitions and their flux 
F int (K km s -1 ) integrated along the whole profile: all 
but H 2 34 S(2 2 o-2n) have been detected. The maps of CepA 
in H 2 S(lio-loi, 2 20 -2n) and SO 2 (3i 3 -2 02 , 5 24 -4 i3 , I6214- 
15 3 i 3 ) had been already reported in a previous paper 
(Codella et al. 2003). The results indicate the occurrence of 
a rich chemistry associated with CepA-East, confirming the 



findings of Bergin et al. (1997). The gas surrounding this re- 
gion is particularly rich in the S-bearing molecules (SO, H 2 S, 
SO2, H 2 CS, HCS+, OCS). Also line spectra due to tracers of 
the high density clumps hosting the YSOs (CS, C ls O, and 
CH 3 C 2 H) as well as to tracers of hot and/or shocked ma- 
terial (SiO, HDO, CH 3 OH) have been obtained. Moreover, 
the use of large bandwidths allowed us to serendipitously 
detect a number of other emission lines due to HC ls O+, 
C 3 H 2 , C 13 CH, CH 2 CO, CH 3 OH, and HNCO, as shown in 
the lower part of Tab. 1. Table 2 reports the lines which have 
not been identified, listing the observed peak frequency, the 
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CepA— East 



CepA East 



o 

0) 

00 

o 
o 



0) 

to 

o 

o 

0) 
Q 




60 30 -30 

R.A. Offset (arcsec) 

Figure 1. Contour maps of the integrated CS J=2-l (upper 
panel) and J=5-4 (lower panel) integrated emission towards 
CepA-East. The maps are overlaid with the grey scale images 
reproducing the integrated emission of H 2 32 S(lio-loi) (upper 
panel) and 32 S02 (3i3-2o2) (lower panel), reported by Codella 
et al. (2003) in their Fig. 1. The empty circles show the IRAM 
beam (HPBW), while the small crosses mark the observed posi- 
tions. The triangles stand for the VLA 2 cm continuum compo- 
nents which trace in the East region two strings of sources arising 
in shocks (Garay et al. 1996). The velocity integration interval 
is —30, +10 km s — 1 . The contours range from 1.0 to 31.0 K km 
s _1 (upper panel) and from 3.0 to 57.0 K km s _1 (lower panel). 
The first contours and the steps correspond to about 3 and 6 a, 
respectively (where a is the r.m.s. of the map). 



spectral resolution, the peak temperature and the FWHM 
linewidth. 

Table 1 shows that the lines detected are associated 
with a wide range of excitation, from fews to hundreds of 
Kelvins. In particular, the detection of emission due to tran- 
sitions above 100 K clearly suggest the presence of high tem- 
perature conditions, which could be associated to the pres- 
ence of a hot core and/or the occurrence of shocked material, 
and consequently requires the analysis of the maps and of 
the line profiles. 
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Figure 2. CS and C S line profiles observed towards CepA East: 
transitions are reported. The dashed line stands for the ambient 
LSR velocity (-10.65 km s" 1 ), according to the CS J = 5-4 and 
C ls O J = 2-1 measurement (see text). 



CepA East 




' LSR 

Figure 3. CepA East: comparison between same line transitions 
of different isotopomers. The solid lines are for the H2 32 S(lio-loi) 
and 32 SO 2 (3i3-2 2) profiles reported by Codella et al. (2003). The 
dotted lines show the corresponding H2 34 S and 34 S02 spectra, 
rescaled to be comparable with the 32 S- lines. Note that only for 
32 SC>2, the 3i3-2o2 line is blended with the I6214-I5313 emission. 
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Table 2. Unidentified detections observed towards CepA-East 



Obs. frequency 


du 




FWHM 


(MHz) 


(MHz) 


(mK) 


(km/s) 


U-102062.23(0.15) 


0.31 


15(3) 


7.8(2.1) 


U-202811. 00(0.19) 


1.00 


96(20) 


5.3(0.9) 


11-213303.65(0.14) 


1.00 


101(12) 


8.4(1.6) 


U-216758.40(0.47) 


1.00 


80(19) 


9.5(1.8) 



3 A VIEW OF CEPA-EAST: DENSE CLUMPS 
AND OUTFLOWS 

3.1 Ambient emission 

We produced relatively small CS maps of CepA-East in 
order to carefully study the distribution of the ambient 
clumps. In particular, we focused our attention on the re- 
gions where the two strings of VLA continuum sources asso- 
ciated with shocks (Garay et al. 1996) and the H2S and SO2 
clumps (Codella et al. 2003) are located. Figure 1 shows the 
maps of the integrated emission due to CS J=2-l (upper 
panel) and J— 5-4 (lower panel) and compare such distri- 
butions with those of H2S Jk_k + = lio-loi (upper panel) 
and SO2 Jk_k + = 3i3-2o2 (lower panel), reported in gray 
scale. The triangles stand for the VLA objects. It is possible 
to note that the two CS maps peak at different positions, 
probably due to the self-absorption of the 2-1 line (see the 
self-reversed profile in Fig. 2). Also the H2S and SO2 emis- 
sions show a different distribution with respect to CS, indi- 
cating that such molecules are not simply reproducing the 
same gas distribution and are able to trace different gas as- 
sociated with different physical and/or chemical conditions. 

Figure 2 reports the CS (upper panels) and C 34 S (lower 
panels) spectra due to the J=2-l and 5-4 transitions. Al- 
though the CS(5-4) line shows the occurrence of wings which 
modify the gaussian profile at about 15% of the line peak, 
such emission can be used to define the characteristics of 
the ambient emission: the LSR velocity is -10.65 km s _1 
(dashed line in Fig. 2), while the FWHM is 3.8 km s -1 . 
The comparison between these CS spectra confirms what 
found in the recent paper by Bottinelli & Williams (2004) 
which clearly indicates that the CS(2-1) profile is affected 
by self-absorption, whereas the CS(5-4) and the two transi- 
tions of C 34 S, due to higher excitation and/or smaller abun- 
dance, are definitely optically thinner and thus well centered 
at the LSR ambient velocity. Moreover, as reported also 
by Bottinelli & Williams (2004), the CS(2-1) line shows a 
blueshifted peak brighter than the redshifted one, suggest- 
ing that infall motions can play an important role in the 
dynamics of the material traced by CS. 

Figure 3 compares the H2 34 S(lio-loi) and the 
34 S02(3i3-2o2) spectra (dotted lines, see also Fig. 1) with 
the profiles of the main isotopomers (continuous lines), re- 
ported in a previous paper (Codella et al. 2003). Note that 
the S02(3i3-2o2) profile is blended with the I6214-I5313 line. 
The comparison between the H2S spectra indicates that the 
double-peak H~2 32 S pattern is caused by self-absorption, with 
the H 2 34 S line peaking at the LSR velocity, where H 2 32 S 
shows the relative minimum. On the other hand, Fig. 3 
confirms that the detected SO2 lines do not present self- 
absorption effects and they peak at the ambient velocity. 



3.2 Kinematics 

Figure 1 shows an elongated structure oriented in the NE- 
SW direction in the CS(5-4) line, which is not affected by 
self-absorption and it has been observed with an angular res- 
olution of 10". To further study this distribution, a CS(5- 
4) channel map is reported in Fig. 4. Focusing the atten- 
tion on the emission at velocities close to the ambient one, 
which corresponds to -10.7 with a FWHM of 1.8 km s~\ 
we note that the elongated structure is due to the presence 
of two clumps centered respectively at the (+20", 0") and 
(-10", -10") offset map positions and with a beam decon- 
volved size of about 15" (0.05 pc at a distance of 725 pc, 
Sargent 1977) at slightly different velocities. This picture 
resembles the distribution of another high-density tracers 
such as N2H + (Bergin et al. 1997): the two clumps could 
indicate star forming sites or could point to gas components 
accelerated by the interaction with the mass loss processes 
indicated by the NE-SW VLA string. The latter case seems 
to be more likely since the clumps are symmetrically located 
with respect to the positions of the driving YSOs (the HW2 
object is located at the centre of the map) and lie along the 
direction corresponding to the main axis of an extremely 
high- velocity outflow («-«lsr>20 km s _1 ; e.g. Narayanan 
& Walker 1996). 

On the other hand, from Fig. 4 it is possible to see 
(i) redshifted (up to -2 km s _1 ) and blueshifted (up to - 
15 km s _1 ) emission associated with the SW jet, and (ii) a 
blueshifted clump located at the end of the SE VLA chain as 
clearly reported in the -15 km s _1 panel. The latter clump 
is also seen in H2S at high velocities (see Fig. 4 of Codella et 
al. 2003) and visible in Fig. 1 (upper panel) at (+40" -20"), 
confirming that at this position we are tracing high density 
material associated with outflow motions. In conclusion, the 
CS emission traces two outflow directions, the SW and SE 
ones, associated with two strings of shocked sources. 

Figure 5 shows the contour maps of the Jk_k + = lio - 
ln, and 2n-2i2 integrated HDO emission towards CepA- 
East. In particular, considering the CS results, the southern 
region has also been mapped to look for signposts of outflow 
motions. The lio-fn emission clearly indicates (i) a struc- 
ture elongated along the SW direction with a peak near the 
centre of the map, and (ii) another elongated feature point- 
ing towards South (see Sect. 3.3). On the other hand, the 
2n-2i2 HDO map, obtained with higher resolution (10"), 
reveals a ~15" (0.05 pc) clump centered at the HW2 coordi- 
nates, where the YSOs are located, and a second unresolved 
structure lying along the SW VLA jet. 

In order to study the HDO kinematics, we present in 
Figures 6 and 7 the velocity channel maps of the I10-I11 and 
2n-2i2 emissions, respectively. From these maps we argue 
that HDO in CepA-East traces the outflow activity as well 
as a small region around the driving YSOs. A quite complex 
scenario is therefore drawn: (i) a central component, clearly 
shown by the 2n-2i2 line, which is associated both with 
ambient velocity and higher redshifted (up to -2 km s _1 ) 
velocities; (ii) a component associated with the SW VLA 
chain and emitting in both HDO lines and (iii) a final com- 
ponent pointing towards South and detected through the 
I10-I11 line. These last two components emit at redshifted 
velocities (up to -4 km s _1 ) and at -12 and -11 km s _1 , i.e. 
at velocities close to the ambient one. Note that the -12 and 
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from 0.14 to 0.56 K km s~ 1 (upper panel) and from 0.60 to 3.00 K 
km s _1 (lower panel). The first contours and the steps correspond 
to about 3 a. 



-11 km s _1 panels of Fig. 6 seem to indicate also the possible 
occurrence of a HDO lio-ln clump at (0",-80"), along the 
southern direction, suggesting a connection with the central 
feature pointing towards South. This southern emission has 
not been detected with the line and therefore for the 

sake of clarity in the corresponding channel maps in Fig. 7 
only a zoom of a smaller region is shown. 

3.3 HDO and H 2 S 

In order to further investigate the southern molecular out- 
flow, the present HDO maps can be compared with the H2S 
ones reported by Codella et al. (2003). In particular, Fig. 8 
shows the contour maps (thin line) of the H2S(lio-loi) emis- 
sion integrated over the redshifted velocity interval typical 
of the southern outflow: -7,-3 km s~ . The HDO(lio-loi) 
map, integrated over the same velocity range, is drawn by 
using a thick contour. From the H2S contours it is possi- 
ble to detect the SE outflow (ending at the position of the 
clump A) and the southern outflow, which extends at the 
position of the clump C. The HDO map is in agreement 
with the picture given by H2S, both suggesting an outflow 
activity towards South. This scenario could reflect different 
chemical and/or physical conditions for the regions where 
the different flows are moving. The southern outflow seems 
chemically enriched with respect to the typical gas compo- 
sition of the dark clouds. It shows a redshifted component 
(see Fig. 6), whereas the blueshifted counterpart could be 
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tentatively singled out in the emission around ~ -12 km 
s~\ i.e. at velocities close to the ambient one, but slightly 
blueshifted with respect to the -10.7 km s _1 value. In partic- 
ular, the tentative detection of HDO emission at (0",-80") 
shown in the channel maps in Fig. 6 is supported by the oc- 
currence of a H2S red- and blueshifted (see Fig. 4 of Codella 
et al. 2003) clump at the same position. 

In Fig. 9 we summarise all the information given by the 
CS and HDO channel maps, drawing a schematic picture 
with the main clumps and outflow directions. The high ex- 
citation HDO(2n-2i2) line is tracing the region close to the 
YSOs, marked in Fig. 9 by an asterisk at (0",0") offset repre- 
senting the HW2 coordinates. Four main flows are identified: 
three (SW, NE, and SE) associated with the VLA sources 
and probably accelerating high density material traced by 
CS. The SW outflow is traced also by HDO emission. In ad- 
dition HDO reveals a fourth outflow pointing towards South 
and not traced by the VLA sources nor by the CS structures. 
Finally, we note that the highest HDO excitation conditions 
occur near the YSOs, where the beam deconvolved line in- 
tensity ratio .Rhdo — T mb (2ii-2i2)/T mb (lio-loi) is aDOut 2. 



On the other hand, the southern flow suggests lower excita- 
tion conditions, with Rnoo — 0.6 at the (0" -20") offset. 



4 A SPECTRAL SURVEY TOWARDS 
CEPA-EAST 

Figure 10 presents the most representative examples of line 
profiles observed towards CepA-East. Different molecules 
exhibit different spectral behaviours and they can be 
grouped in three classes, summarised in Table 3: (i) am- 
bient species (C 13 CH, C 3 H 2 , CH 2 CO, CH3C2H, HC ls O+), 
not detected at high velocities and with relatively narrow 
(3-4 km s" 1 ) lines at velocities close to the ambient one 
(hereafter called ambient peak), (ii) outflow tracers (CO, 
CS, H 2 S, SiO, SO, SO2), which show extended wings (e.g. 
v-fLSR up to ~ 30 km s _1 for CO and 10 km s _1 for SO) and 
span the whole range of observed velocities, and (iii) species 
(OCS, H 2 CS, HDO, and CH3OH) which are associated with 
wings and, in addition, show a redshifted secondary peak at 
-5.5 km s _1 (hereafter called outflow peak), well separated 
(by about 5 km s" 1 ) from the ambient velocity. Unfortu- 
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Figure 7. Zoom of the channel map of the HDO Jk_k + = 2n— 2i2 emission towards CepA-East. Each panel shows the emission 
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Table 3. Summary of the properties of the profiles observed towards CepA-East 



Molecular tracers in CepA-East 

CH 3 OH, HCS+, HDO, H 2 CS, OCS =>• Ambient, Outflow, and Peak @ -5.5 km s" 1 

CO, CS, H 2 S, SiO, SO, SO2 => Ambient and Outflow 

C 13 CH, C3H2, CH2CO, CH3C2H, HC ls O+ =^ Only Ambient 



nately, the present data do not allow us to clarify the spatial 
distribution of the spectral outflow peak, i.e. whether it is 
tracing a small clump or it is related to a more extended 
structure. However, the comparison between HDO and H2S 
emissions shown in Fig. 9 clearly suggests that the spectral 
outflow peak is tracing an elongated structure flowing to- 
wards South. Only accurate CH3OH, H2CS, and OCS maps 
can confirm this conclusion. 

Figure 11 displays the CH 3 OH(5 K -4 K ) and (3 k -2k) 
spectra in the upper and lower panel, respectively: the 
dashed lines indicate the predicted positions of the hyper- 
fine component at different excitations. Although different 
lines are blended, from these spectacular emissions it is pos- 
sible to see that almost all the components, including those 
at very high excitations with E u ~ 120 K, have been de- 
tected and show the redshifted peak. It is worth noting that 
other very high excitation CH3OH lines (up to E u = 390 
K) have been serendipitously detected, as reported in Table 
1. Finally, Fig. 12 shows the striking differences observed 
on the profiles of several transitions of HDO (upper pan- 



els) and H 2 CS (lower panels). From these spectra, it is clear 
that the two HDO peaks are associated with roughly sim- 
ilar excitation conditions, whereas for H2CS, the -5.5 km 
s~ peak increases its intensity with respect to the -10.7 km 
s _1 peak with excitation indicating different physical condi- 
tions. For instance, the ratio between the mean brightness 
temperatures of the outflow and ambient peaks is ~ 0.6 for 
the Jk_k + = 7i6-6i5 transition at 60 K, while it is < 0.1 
for the Jk_k + = 3i3-2i2 line at 23 K. 

In conclusion, the obtained line profiles indicate that 
at the CepA-East positions the emissions due to the clump 
hosting the YSOs and to the associated molecular outflows 
coexist and produce two distinct line peaks which suggest 
different excitation conditions. The results suggest also that 
different molecules can trace different excitation conditions 
at the same observed velocity. This gives us a precious op- 
portunity to perform a multiline analysis in order to clarify 
the physical conditions associated with the ambient emis- 
sion and those associated with the southern chemically rich 
molecular outflow. 
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Figure 8. Contour maps (thin line) of the H2S(lio-loi) red- 
shifted emission (Codella et al. 2003) integrated over the velocity 
interval —7,-3 km s _1 . Two main flows are clearly drawn: one 
pointing to SE (and ending at the clump A position), and an- 
other towards South (ending at clump C). Symbols arc drawn 
as in Fig. 1. The thick solid line shows the contour maps of the 
HDO(lio-loi) emission integrated over the same range of veloc- 
ity. The thick dotted line shows the distribution of the HDO(lio- 
loi) emission integrated over the interval -13,-11 km s" 1 and 
tentatively points out a blueshifted counterpart (see the channel 
maps of Fig. 6). 



5 DERIVED GAS PARAMETERS 

5.1 Analysis procedures 

By means of statistical-equilibrium calculations and using 
the four observed SO lines, it is possible to estimate the total 
column density (N to t) as well as the kinetic temperature (Tk) 
and the hydrogen density (riH 2 ). A Large Velocity Gradient 
(LVG) model and the collisional rates from Green (1994) 
have been used. Due to the lack of SO maps, these calcula- 
tions have been performed using the mean-beam brightness 
temperatures relative to the intensities of the four spectra, 
without beam filling factor correction. If the source was def- 
initely smaller than the four beamwidths, the LVG calcula- 
tions would lead to an overestimate of the column densities 
(up to a factor 9) and of excitation conditions (up to factors 

1.2 and 2 for Tk and bh 2 , respectively). 

Following Cesaroni et al. (1991), statistical equilibrium 
computations in the LVG approximation have been used to 
analyse the CS and C 34 S spectra. In particular, for opti- 
cally thin conditions, the line brightness temperatures de- 
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Figure 9. A schematic picture (not to scale) of the directions of 
the multiple outflows driven by the CepA-East YSOs as traced by 
HDO and CS emission (see text). Following the channel maps four 
main flows are identified: three (SW, NE, and SE) associated with 
the VLA shocked chains and with CS clumps, and one (S) traced 
by the HDO(lin-loi) emission. The higher excitation HDO(2n- 
2i2) line well defines the SW outflow and with the central region 
close to the YSOs coordinates (marked by an asterisk). 



pend only on density and, to a minor degree, on tempera- 
ture. In this case, we excluded CS(2-1) which clearly shows 
self-absorption effects and we used CS(5-4) as well as the 
C 34 S lines to derive 71h 2 estimates. 

For the CH 3 OH, CH3C2H, H 2 CS, and 34 S0 2 molecules, 
observed through at least three lines, in order to estimate the 
rotational temperature and column densities, the standard 
rotation diagram method, assuming LTE and optically thin 
conditions, has been used. For the other molecular species, 
the total column densities have been calculated using the 
constants given in the databases for molecular spectroscopy 
and considering the temperature estimates derived from the 
LVG and rotation diagram results. For comparison, also 
the values calculated by using the H2S and SO2 emissions 
(Codella et al. 2003) are reported. 

When possible, different column densities have been de- 
rived for the different lines at -10.7 and -5.5 km s _1 . In 
particular, the methanol patterns clearly allow us to derive 
different excitation conditions, while H2CS leads to an esti- 
mate of the temperature for the -10.7 km s" 1 component. 
Finally, following the SO2 analysis of Codella et al. (2003), 
we have assumed two components also for CS and SO which 
show strong redshifted wing emission with an intense emis- 
sion around -5 km s" 1 . 



5.2 Results 

The derived parameters have been summarised in Table 4. 
The LVG results based on the SO spectra lead to kinetic 
temperatures of 60-100 K and iih 2 ~ 5 10 6 cm" 3 for the 
-10.7 km s _1 line at ambient velocity, while slightly higher 
temperatures (70-180 K) and definitely higher densities (2 
10 6 -6 10 7 cm ) have been derived for the -5.5 km s com- 
ponent. The total SO column densities are around 4 10 14 
and 5 10 13 cm" 2 for the ambient and outflow components, 
respectively. 
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Figure 10. Molecular line profiles observed towards CepA East: 
species and transitions are reported. The dashed lines stand for 
the ambient LSR velocity (-10.65 km s^ 1 ), according to the C ls O 
measurement, and for the component at -5.51 km s _1 well out- 
lined e.g. in the H2CS profile. In case of the CH30H(5k — 4k) 
spectrum, four hyperfine lines at different excitations are present: 
K=0 A+ (E u =35 K), K=-l E (£ u =39 K), K=-4 E (£ u =122 K), 
and K=A A+" (B u =115 K). The CH 3 C 2 H(6k-5k) pattern shows 
the #=0,1,2,3 lines. The small vertical labels CH3OH (£ u =126 
K), C3H2, and HNCO refers to three serendipity detections (see 
Table 1). 



The CS LVG results indicate hydrogen number densi- 
ties larger than 6 10 4 cm~ 3 for the ambient component. For 
the outflow component, the CS analysis leads to densities in 
the range between 4 10 3 and 5 10 4 cm -3 . On the other hand, 
the C 34 S data for the ambient component suggest lower den- 
sities with values around 2 10 4 cm~ 3 . The C 34 S emission is 
optically thinner than the CS one, and in fact CS shows self- 
absorption whereas C 34 S shows Gaussian profiles; thus, it is 
reasonable to expect that C 34 S is tracing inner regions with 
respect to CS. In this case, the LVG results seem to give puz- 
zling results for CepA-East with a high density envelope and 
a lower density central region. One possible solution comes 
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Figure 11. CH 3 OH(5 K -4 K ) and CH 3 OH(3 K -2 K ) line profiles 
observed towards CepA East. The dotted lines stand for the 
predicted positions of the hyperfine components at different ex- 
citations. The HNCO(llon-lOoio) line is also marked in the 
methanol (5k-4k) pattern. The velocity scale of the spectra is cal- 
culated with respect to the frequency of the 5o-4rjA + (241791.43 
MHz) and 3 -2 A+ (145103.23 MHz) lines. 



from the recent results reported by Bottinelli & Williams 
(2004), who analyse the large scale dynamics of CepA-East 
by using a density profile measured from a 850 fim map of 
the region. The model proposed by these authors has a high 
density center with depleted CS and an outer envelope with 
a low nu 2 and normal CS abundance. Therefore, our hydro- 
gen density estimates could be biased by the fact that the 
CS abundance is varying along the central region. 

The least-square fits to all detected lines in the CH3OH 
rotation diagram, shown in Fig. 13, gives, for the ambient 
component (middle panel), rotational temperatures of about 
30 K, while the total methanol column density is about 7 
10 14 cm~ 2 . The CH3OH temperature is confirmed by the 
CH3C2H rotation diagram which leads to 37 K. This sug- 
gests that at the ambient velocities both molecules trace 
similar environments. On the other hand, for the outflow 
component (lower panel of Fig. 13) the derived tempera- 
ture is 220 K and the CH3OH column density is 3 10 15 
cm~ 2 . However, due to the blending between ambient and 
outflow components for the methanol spectral patterns, the 
plot for the outflow component is less clear and the fit is 
definitely more doubtful. In any case, the comparison be- 
tween the plots of the ambient component and the outflow 
one, which shows high excitation (E u = 390 K) CH3OH 
emission, clearly indicates a definitely higher temperature 
for the outflow component. 

Using the H2CS and 34 S02 emissions, an estimate of the 
rotational temperature of the ambient emission, observed in 
three lines, has also been obtained. The results lead to low 
temperatures (20 K) for H2CS and definitely higher values 
(200 K) for 34 S0 2 , suggesting that the H 2 CS spectral peak 
at rest velocity mainly traces relatively cool ambient mate- 
rial. The column densities are Nn 2 cs — 3 10 13 and N3i S o 2 
- li) :; cm '. 

Taking into account the temperature estimates obtained 
from the SO and CH3OH data, the total column densities for 
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Table 4. Column densities, temperatures, and hydrogen densities for the ambient and outflow gas components 
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a The column densities and the abundances regarding the species detected through only one line have been derived assuming 
X[CO]/X[H 2 ] = 10- 4 , and T^\ n = 20-100 K (ambient peak) and T^\ n = 200 K (outflow peak), following the temperature measurements 

based on the CH3C 2 H, CH3OH, SO, and S0 2 emissions. b Average ratio between the outflow and the ambient abundances. c For 
C 13 CH and HC 18 0+, the abundances refer to the most abundant isotopomcrs (CCH and HCO+) calculated by assuming 12 C/ 13 C = 
89 and 16 0/ 18 = 490 (Wilson & Rood 1994). d For H 2 S and SiO the ambient and outflow components cannot be disentangled; thus 
the values refer to the whole emission and the abundances arc not reported. e The H 2 S and S0 2 results are taken from the Codella et 
al. (2003) paper: the H 2 34 S emission has been used to derive the optical depth and consequently refine the H 2 S column density. 



the molecular species observed in one or two lines have been 
calculated using the spectroscopic constants given in the lit- 
erature and summarised in the JPL (Pickett et al. 1998), 
NIST (Lovas 2004), Cologne (Miiller et al. 2001), and Lei- 
den (Schoier et al. 2005) databases and considering a kinetic 
temperature in the range 20-100 K for the ambient compo- 
nent and 200 K for the outflow one. Table 4 reports the 
values of the total column densities for the ambient emis- 
sion: A^ c i3ch an d Nc 3 n 2 are about 1 10 14 cm -2 , A^hnco 
and A?ocs ^ 3 10 13 cm" 2 , N C n 2 co ~ 5 10 15 cm" 2 , N HG i So+ 
and iV HCS+ ~ 3-5 10 12 cm" 2 , while Nudo is around 3 10 14 
cm- 2 . On the other hand, for the outflow component we 
have column densities around 3 10 12 , 6 10 14 , and 5 10 13 
cm -2 for HCS+, HDO, and OCS, respectively. We note that 
the HCS + /CS abundance ratio is quite large, ~ 0.1, in the 
outflow peak, close to that measured in dark clouds (Ohishi 
et al. 1992), suggesting that the HCS + emission has been 
overestimated, due the low S/N ratio of the HCS + outflow 
peak, and/or that the derived CS abundance is underesti- 
mated, probably due to the blending between ambient and 
outflow spectral components (Sect. 5.1). 

The SiO results shown in Table 4 refer to the multiline 
SiO survey given by Codella et al. (1999), whereas the H 2 S 
and S0 2 results are taken from the Codella et al. (2003) 
paper: the H 2 34 S emission have been used here to derive 
the optical depth, T(i 10 -i 01 ) — 1-2, and consequently refine 
the H 2 S column density. Finally, it is worth noting that for 
H 2 S and SiO the ambient and outflow components cannot 
be disentangled both showing continuous extended wings, 
and thus the values in Table 4 refer to the whole emission. 



In conclusion, we find that in CepA-East at ambient 
velocities (i) the gas is associated with high densities, > 
10 5 cm -3 , and (ii) different components at different tem- 
peratures coexist, ranging from the relatively low kinetic 
temperatures, less than 50 K, measured with H 2 S, CH3OH, 
H 2 CS, and CH3C2H, to definitely higher temperature con- 
ditions, ~ 100-200 K, obtained from the SiO, SO, and S0 2 
spectra. In other words, at ambient velocities we are prob- 
ably sampling different layers of the high-density medium 
hosting the star forming process, some of them probably 
heated by the stellar radiation. On the other hand, for the 
outflow component we derive two density regimes: about 10 4 
cm- 3 from CS and around 10 7 cm- 3 from SO, whereas the 
temperatures are always quite high: ~ 100-200 K, indicating 
regions compressed and heated by shocks. 



6 DISCUSSION 
6.1 Abundances 

With the aim of obtaining an estimate of the abundances 
of the observed species, the H 2 column densities have been 
derived from C 1 O for the ambient component and from 
CO and C ls O for the outflow peak, by assuming a standard 
X[CO]/X[R 2 ] ratio (~ 10" 4 ; Frerking et al. 1982, Lacy et al. 
1994). We obtained the fractional abundances with respect 
to H 2 as the ratio between column densities and finally an 
enhancement factor / for each molecule has been derived as 
the ratio between the abundances in the outflow and in the 
ambient peaks. The derived values are reported in the last 
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Figure 12. HDO (upper panels) and H2CS (lower panels) line 
profiles observed towards CepA East: transitions arc reported (see 
Table 1). Note how for H2CS the -5.5 km s _1 component in- 
creases its intensity with respect the -10.7 km s peak with 
excitation. 



column of Table 4, where two classes can be distinguished. 
On the one hand, CS, SO, and SO2 do not show a definite 
abundance enhancement (/ ~ 1), while the HCS + abun- 
dance increases only by a factor of 4. On the other hand, 
CH3OH, HDO, H 2 CS, and OCS appear to be definitely en- 
hanced by at least one order of magnitude (/ ~ 10-20), 
confirming that these molecules are closely associated with 
the chemistry occurring in molecular outflows (Bachiller & 
Perez Gutierrez 1999). 

Given the well known observational difficulties of de- 
tecting the H2O isotopomer, HDO represents with H3, 8 an 
alternative tool to investigate the occurrence of water emis- 
sion in space. In fact, HDO has been already observed to- 
wards high-mass star forming regions, mainly hot cores (e.g. 
Jacq et al. 1990, Helmich et al. 1996, Gensheimer et al. 1996 
and references therein). In addition, the high velocity peak 
indicates that HDO, and thus water, is a very sensitive shock 
tracer in molecular outflows driven by YSOs. The present 
data do not allow us to give a direct estimate of the H 2 
abundance. However, by assuming the cosmic [D]/[H] ra- 
tio, ~ 1.5 10~ 5 (Oliveira et al. 2003), or the typical [D]/[H] 
ratios derived for hot cores, 2-6 10~ 4 (Jacq et al. 1990), 
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Figure 13. Rotation diagrams for the CH3C2H (upper panel) 
and CH3OH (middle and lower panels) transitions measured to- 
wards CepA-East. For the methanol emission two spectral regimes 
have been considered: ambient and outflow (see text). The pa- 
rameters 7V U , g u , and E u are, respectively, the column density, 
the degeneracy and the energy for the upper levels of the transi- 
tions reported. The derived values of the rotational temperature 
are reported. 



it is possible to roughly calculate a H 2 abundance from 
our HDO measurements. We therefore derive for CepA-East 
X[H 2 0]/X[H 2 ] in the 10~ 6 -10~ 4 range for the ambient peak 
and ~ 10~ 5 -10~ 4 for the outflow peak, confirming that the 
water abundance can be extremely enhanced in star forming 
regions and that H2O is a major coolant in warm gas (van 
Dishoeck & Blake 1998 and references therein). 

Finally, the high abundances of CH3OH, HDO, H 2 CS, 
and OCS, which show a well separated high velocity range 
(see Figs. 10 and 12), are useful tools to investigate the 
shocked material by comparing to those of other shock trac- 
ers such as SiO, SO, and SO2. 



6.2 Different tracers at different velocities 

In order to investigate the properties of the line wings which 
do not show the peak at -5.5 km s~ , Fig. 14 reports the 
distribution with velocity of the brightness temperature ra- 
tio between SO(7-6) (upper panel), S02(524~4i3) (middle 
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Figure 14. Distribution with velocity of the ratio between the 
brightness temperatures of the SO(7-6) (E u = 35 K; upper panel), 
H 2 S(lio-loi) {Eu = 28 K; middle panel), S0 2 (5 2 4-4i 3 ) (E u = 
24 K; lower panel) lines and of the SiO(5-4) {E u = 21 K) profile 
as detected towards CepA-East. The H 2 S and S0 2 data are taken 
from Codclla ct al. (2003). The HPBWs of the observations used 
here are: 9" (SO), 10" (S0 2 ), and 11" (H 2 S, SiO). The dot- 
dashed line underlines the ambient LSR velocity according to the 
CS J = 5-4 and C ls O J = 2-1 measurement. 

panel), H 2 S(lio-loi) (lower panel) and SiO(5-4), as ob- 
served towards CepA-East. Note that for such a comparison 
we selected transitions observed with a similar HPBW (SO: 
9", S0 2 : 10", H 2 S and SiO: 11") and with similar excitation 
{E u ~ 21-35 K, see Table 1). From Fig. 14 it is possible to 
see that: 

(i) SiO emission dominates at the largest velocities, where 
the highest excitation conditions are expected. Actually, a 
comparison between the line profiles of different excitation 
transitions of H 2 S and S0 2 observed towards CepA-East 
clearly indicates that the higher the velocity the higher the 
excitation (see Fig. 8 of Codella et al. 2003) ; 

(ii) among the S-bearing species, S0 2 has the distribution 
associated with the largest width, followed by SO, and finally 
by H 2 S. This is confirmed also by the SO/H 2 S, S0 2 /H 2 S, 
and S0 2 /SO column density ratio profiles, not shown here, 
which increase with velocity. 

These results confirm the close association of SiO with 
the shocks occurring along the molecular outflows. Sputter- 
ing on dust grains is probably the most efficient mechanism 
that leads to a SiO abundance enhancement, although also 
grain-grain collisions can play a role depending on the gas 
density (Schilke et al. 1997, Caselli et al. 1997). However, 
it seems that H 2 S, SO, and S0 2 preferentially trace more 
quiescent regions than SiO. from the different widths of the 
resultant column density ratio profiles, we note the lack of 
an intense H 2 S emission at the highest velocities. A possi- 
ble explanation is a fast conversion of H 2 S into SO and S0 2 



and possibly OCS, in agreement with the chemical models 
where first H 2 S is injected from grains, then the other S- 
bearing species are formed very quickly in ~10 3 yr (Pineau 
des Forets et al. 1993, Charnley 1997, Viti et al. 2004). Alter- 
natively, it may be that hydrogen sulphide is not the major 
sulphur carrier in the grain mantles, as suggested also by the 
lack of H 2 S features in the ISO spectra (Gibb et al. 2000, 
Boogert et al. 2000) and by the low temperature derived 
from the present H 2 S spectra, ~ 27 K, well in agreement 
with that measured by van der Tak et al. (2003, 25 K) to- 
wards high-mass star forming regions. 

On the other hand, the occurrence of the outflow peak 
at -5.5 km s _1 shows that: 

(i) OCS and H 2 CS emit at high velocities, where SiO 
emission starts to dominate and the excitation conditions 
are high. At these velocities also the abundance of CH3OH, 
which is expected to be abundant in grain mantles and to be 
released in the gas phase after shocks, definitely increases. 
These results suggest that OCS and H 2 CS are among the 
sulphur carriers in the grain mantles, or that at least they 
are rapidly formed once the mantle is evaporated. Both sce- 
narios are in agreement with the results by van der Tak et al. 
(2003) who measured high excitation temperatures for OCS 
(~ 100 K). Unfortunately, the present data do not allow us 
to derive a direct temperature estimate from the OCS data. 
However, an estimate of the temperature of the gas traced 
by OCS is given by the CH3OH emission, which well defines 
the OCS outflow peak and leads to a high temperature: 220 
K; 

(ii) HDO shows an enhancement of its fractional abun- 
dance (X[HDO]/X[H 2 ]) at such high velocities. From the 
estimated value, 10 -8 , we can roughly infer the abundance 
of water assuming a D/H ratio equal to 1.5 10~ 5 (Oliveira 
ct al. 2003): X[H 2 0]/A:[H 2 ] ~ 6 10~ 4 . This suggests that 
the majority of oxygen is here locked into water molecules 
(Meyer et al. 1998). This is expected if we are observing 
regions heated by shocks above 200 K, so that the major- 
ity of OH molecules quickly produce water through reac- 
tions with molecular hydrogen (e.g. Hartquist et al. 1980). 
The methanol measurements confirm the occurrence of high 
temperatures at these velocities. 

Finally, note that sulphur could be also in atomic form 
when evaporated from the grain mantles, as suggested by 
Wakelam et al. (2004). However, the atomic S is expected 
to be very quickly (< 10 3 yr) locked into SO and S0 2 , so 
that a big amount of atomic S is not expected to last in 
the gas phase for a long time and thus difficult to observe. 
The form that sulphur takes on dust grains is still far to 
be known: the present data indicate that OCS and H 2 CS 
may play an important role and/or that they are effectively 
formed in high-temperature gas. 

In conclusion, the present analysis of the line profiles 
indicates that H 2 S, SO, and S0 2 may not be easily used 
as chemical clocks of the shocked material. On the other 
hand, H 2 CS and OCS could be useful candidates: in the 
next section, a possible use of their emission is discussed. 

6.3 On the origin of the outflow spectral peak 

Since the outflow peak is clearly defined by two sulphuret- 
ted species like OCS and H 2 CS and by two shock tracers 
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Figure 15. Evolution of the abundance ratios H2CS/OCS (upper panels) and CH3OH/H2CS (lower panels) as functions of time for a 
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like CH3OH and HDO, we attempt to use such emission as 
a chemical clock. We compared the observations with the 
theoretical calculations recently reported by Wakelam et al. 
(2004). The authors developed a time-dependent chemical 
model of sulphur chemistry in hot-cores with up-to-date 
reaction rate coefficients, following the molecular compo- 
sition after the injection of grain mantle species into the 
gas phase. They found that the abundances of the main 
S-bearing species (H2S, SO, OCS, and SO2) strongly de- 
pend on the physical conditions, on the oxygen abundance 
in the gas phase, on the adopted grain mantle composition 
as well as on the time. Hence, the use of the abundance 
ratios can be often useless to derive age estimates. Never- 
theless, Wakelam et al. (2004) compared the observed and 
predicted abundance ratios for the hot cores Orion KL and 
IRAS16293-2422, for which the physical conditions had been 
previously obtained. The authors were able to reproduce the 
observed abundance and give an age estimate only assuming 
that a large amount of atomic sulphur is initially present in 
the post-evaporative gas, with X[S]/X[H2] between 3 10~ 5 
and 3 10~ 6 . 

In view of these results, we applied the Wakelam et al. 
(2004) model to the observed emission of the CepA out- 
flow peak. Although the model has been initially used to 
investigate hot cores, it follows how the S-bearing molecu- 
lar abundances vary with time when the gas undergoes a 
sudden change in its temperature and density, and in its 



overall chemical abundance, because of the evaporation of 
grain mantles. Therefore, in principle, the model can be rea- 
sonably used also for shocked conditions occurring along 
a molecular outflow. However, the physical conditions and 
the evolution with time should be different from the hot 
core case since the physical processes are different. Magneto- 
Hydrodynamics models predict that during the passage of a 
shock wave, the gas temperature and density simultaneously 
increase before decreasing after the shock passage on time 
scales which depend on the considered model (e.g. Bergin et 
al. 1998, Flower & Pineau des Forets 2003). In our cases, the 
physical conditions are difficult to constrain in the shocked 
regions because we do not know how the physical conditions 
evolved and also because the observed molecules indicates 
temperature and density ranges (see Table 4). For these rea- 
sons, we considered three different models. In the first two 
models, the temperature and the density suddenly increase 
to 1000 K and 10 7 cm -3 and the chemistry evolves for 100 
(Case 1) and 1000 yr (Case 2), respectively. After the pas- 
sage of the shock, the temperature decreases to the observed 
values, i.e. ~ 200 K. The chemistry of the post-shocked gas 
then evolves for 10 6 yr. In the third model (Case 3), we as- 
sumed that the observed molecules mostly trace the outer 
material of the shock where the temperature and the density 
only increase to the observed values without passing through 
a phase with temperatures around 1000 K. This last scenario 
is motivated by the fact that at the CepA distance (725 pc), 
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the filling factor relative to the shocked regions where the 
temperature is expected to exceed 1000 K could be definitely 
small and thus their emission could be strongly diluted. 

From a chemical point of view, the starting point is the 
dark molecular cloud composition computed by the chemical 
model (see details of composition A in Wakelam et al. 2004) . 
Due to the shock passage, the molecules contained on the 
grain mantles, such as H2O, H2CO, CH3OH and S-bearing 
species, are sputtered in the gas phase. For the other species, 
we took the abundances observed as icy features towards 
massive star forming regions (see Wakelam et al. 2004). 

Since the initial form and abundance of sulphur sput- 
tered from grains is still an open question, as a first step, we 
decided to use the results of Wakelam et al. (2004) who re- 
produce reasonably well the S-bearing observations towards 
hot cores using a model in which sulphur is mainly evap- 
orated from grains in atomic, OCS, and H2S forms. The 
initial fractional abundances are 10 -7 for OCS and H2S, 3 
10~ 6 for the S form (Palumbo et al. 1997, van Dishoeck & 
Blake 1998, Wakelam et al. 2004). The implicit assumption 
here is that the majority of sulphur is depleted in the refrac- 
tory grain cores. 

Since HDO is not included in the Wakelam et al. (2004) 
model, we only used the abundance ratios between CH3OH, 
OCS, and H2CS. Figure 15 reports the evolution of the 
H2CS/OCS (upper panels) and CH3OH/H2CS (lower pan- 
els) ratios as a function of time for the Case 1 (shocked 
gas 100 yr old). The initial abundances of the three species 
are 9.8 10~ 8 , 8.20 10 -11 , and 3.9 10~ 6 , for OCS, H 2 CS, 
and CH3OH, respectively. Following the physical parame- 
ters derived for the outflow peak (see Table 4) and adopt- 
ing a conservative approach, we used temperatures of 100 
K (left panels) and 300 K (right panels), whereas for the 
hydrogen density we assumed 10 7 cm~ 3 (dot-dashed line). 
The horizontal thick lines indicate the values derived from 
the observations towards CepA-East for the -5.5 km s _1 
peak: 2 for H 2 CS/OCS and 30 for CH3OH/H2CS. Figure 15 
shows that, for the assumed nu 2 and Tki n ranges, the ob- 
served ratios lead to ages for the shocked gas in the CepA- 
East redshifted component in the range ~ 10 3 -3 10 4 yr. The 
H2CS/OCS trend is due to a maximum of the H2CS abun- 
dance at 10 4 yr and to a decrease of the OCS abundance 
after 10 3 yr. On the other hand, the CH3OH/H2CS trend is 
due to a dramatic decrease of the methanol abundance for 
ages > 10 4 yr. Note that the CH3OH/H2CS result solves 
the dichotomy given by the H2CS/OCS ratio, which is less 
selective, proposing ages either < 10 4 yr or ~ 10 6 yr. 

Cases 2 and 3 (shocked gas 1000 yr and warmed dark 
cloud) do not show significant changes and thus lead to the 
same conclusion, leaving open the question about the ori- 
gin of the outflow peak: gas which passed through a hot 
phase or warm gas in the shock surroundings? In any case, 
even taking into account all the uncertainties of the col- 
umn density estimates, of the abundance ratios, and of the 
physical conditions, the comparison between the observa- 
tions and the model used here allow us to obtain an age 
for the shocked gas in CepA-East of the order of ~ 10 3 -3 
10 4 yr. Only precise measurements of density and temper- 
ature, which are usually hampered by the high degree of 
confusion associated with the observed regions, can refine 
age estimates further. Future investigations of emission due 
to S-bearing species towards other molecular outflows and 



hot cores where the physical conditions are known will con- 
firm whether the CH3OH/H2CS and H 2 CS/OCS abundance 
ratios can in fact be used to estimate ages. 



7 SUMMARY 

The CepA star forming region has been investigated through 
a multiline survey at mm-wavelengths. The results indicate 
the occurrence of a rich chemistry surrounding the YSOs of 
the CepA-East stellar association. The main findings are as 
follows: 

1. The CS and HDO maps draw a complex scenario, de- 
tecting high excitation clumps hosting the YSOs driving 
multiple outflows. Four main flows have been identified: 
three are along the SW, NE, and SE directions, traced also 
by strings of VLA continuum sources, and are accelerating 
high density CS clumps. In addition, HDO reveals a fourth 
outflow pointing towards South, which had been previously 
detected only through H2S and SO2 observations, and is 
then associated with conditions particularly favourable to a 
chemical enrichment. 

2. At the CepA-East position different molecules exhibit 
different spectral behaviours: three classes can be identified. 
Some species (C 13 CH, C 3 H 2 , CH 2 CO, CH3C2H, HC 18 0+) 
peak with relatively narrow lines at ambient velocities (am- 
bient peak). Other molecules (CO, CS, H 2 S, SiO, SO, SO2) 
show extended wings and trace the whole range of the out- 
flow velocities. Moreover, there is a group of species (OCS, 
H2CS, HDO, and CH3OH) which shows wings and, in ad- 
dition, well defines a high velocity redshifted spectral peak 
(outflow peak) which can be used to investigate the SE-S 
outflows. 

3. The physical conditions associated with both the am- 
bient and outflow peaks have been estimated. By means 
of statistical-equilibrium calculations, using LVG codes, we 
have obtained the physical parameters of the molecular gas 
traced by SO and CS. For CH3C2H, H 2 CS, 34 S0 2 , and 
CH3OH we have used the rotation diagram method. For 
the other molecular species, we have calculated the total 
column densities assuming LTE conditions. At ambient ve- 
locities the gas is quite dense (> 10 5 cm~ 3 ) and different 
components at different temperatures coexist, ranging from 
the relatively low kinetic temperatures (< 50 K) measured 
with H 2 S, CH3OH, H 2 CS, and CH3C2H, to definitely higher 
temperature conditions, ~ 100-200 K, obtained from the 
SiO, SO, and SO2 spectra, which may trace layers directly 
heated by the stellar radiation. For the outflow component 
densities between ~ 10 4 cm~ 3 to ~ 10 7 cm~ 3 and high tem- 
peratures, ~ 100-200 K, have been found, indicating regions 
compressed and heated by shocks. 

4. The comparison between the line profiles of different 
outflow tracers of molecular outflows shows that SiO domi- 
nates at the highest velocities, where the highest excitation 
conditions are found. This confirms the close association of 
SiO with shocks. On the other hand, H2S, SO2, and SO 
preferentially trace more quiescent regions. In particular, we 
find a lack of a bright H2S emission at the highest veloci- 
ties. Moreover, OCS and H2CS emit at quite high velocities, 
where (i) SiO emission dominates and the excitation con- 
ditions are high, and (ii) CH3OH and HDO, other shock 
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tracers, increase their abundance. These results could indi- 
cate that H2S is not the only major sulphur carrier in the 
grain mantles, and that OCS and H2CS may probably play 
an important role on the grains; or that alternatively they 
could rapidly form once the mantle is evaporated after the 
passage of a shock. 

5. We checked the possible use of the CH3OH, OCS, and 
H2CS emission as chemical clocks to measure the age of 
the shocked material associated with the outflow peak. The 
observations have been compared with the theoretical cal- 
culations recently reported by Wakelam et al. (2004), who 
developed a time-dependent model with up-to-date reaction 
rate coefficients for the sulphur chemistry. Once associated 
with the derived physical parameters, the H2CS/OCS and 
CH3OH/H2CS column density ratios led to ages in the range 
10 3 -3 10 4 yr. Further observations of emission due to sul- 
phuretted molecules as well as accurate measurements of 
density and temperatures associated with the outflowing 
material are necessary in order to calibrate these chemical 
clocks. 
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